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ABSTRACT 

We study the local stability of stratified, differentially-rotating fluids to 
axisymmetric perturbations in the presence of a weak magnetic field and of 



> , 

O , finite resistivity, viscosity and heat conductivity. This is a generalization of the 

Goldreich-Schubert-Fricke (GSF) double-diffusive analysis to the magnetized 
^ ! and resistive, triple-diffusive case. Our fifth-order dispersion relation admits a 

Tj- ! novel branch which describes a magnetized version of multi-diffusive modes. 

^ ' We derive necessary conditions for axisymmetric stability in the inviscid and 



perfect-conductor (double-diffusive) limits. In each case, rotation must be 
constant on cylinders and angular velocity must not decrease with distance from 
the rotation axis for stability, irrespective of the relative strength of viscous. 



o 

! resistive and heat diffusion. Therefore, in both double-diffusive limits, solid body 

rotation marginally satisfies our stability criteria. The role of weak magnetic 
^ i fields is essential to reach these conclusions. The triple-diffusive situation is 

^ ! more complex, and its stability criteria are not easily stated. Numerical analysis 

of our general dispersion relation confirms our analytic double-diffusive criteria, 
but also shows that an unstable double-diffusive situation can be significantly 
stabihzed by the addition of a third, ostensibly weaker, diffusion process. We 
describe a numerical application to the Sun's upper radiative zone and establish 
that it would be subject to unstable mult i- diffusive modes if moderate or strong 
radial gradients of angular velocity were present. 

Subject headings: accretion disks — hydrodynamics — MHD — instabilities — 
turbulence — Sun: rotation, interior, magnetic fields — stars: rotation 
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1. Introduction 

Developments in the last decade have made it clear that magnetic fields, even weak 
magnetic fields, are essential to our understanding of the dynamics of differentially rotating 
accretion disks and accretion flows in general (Balbus & Hawley 1991; 1998; Balbus 2003; 
Blaes 2003). Although Keplerian rotation profiles are hnearly stable to hydrodynamical 
axisymmetric perturbations, the introduction of a weak magnetic field, acting as a tether 
between fiuid elements, renders such disks unstable to the magnetorotational instability, 
or MRI. The ensuing MHD turbulence is now generally viewed as the primary mechanism 
providing the outward angular momentum transport responsible for accretion in sufficiently 
ionized, non self-gravitating disks (Hawley, Gammie & Balbus 1995; Armitage 1998; Hawley 
2001; Balbus & Hawley 1998). 

The destabilizing effect on stellar differential rotation of an embedded magnetic field 
was studied well before the accretion disk community took notice of its importance. Fricke 
(1968) pointed out, for example, that time-steady field configurations (isovelocity rotation 
contours running along field lines) need not be stable. Acheson's thoroughgoing and detailed 
review (1973) analyzed the effects of toroidal fields, noting that even very weak fields 
quahtatively change the Goldreich-Schubert-Pricke (Goldreich & Schubert 1967; Fricke 
1968) criterion for rotational stability, a topic with which we shall be very much concerned 
in this paper. The magnetic field problem in its full generality was avoided, however, 
because of the apparent complications associated with a time-dependent equilibrium field 
caused by shear. The clear understanding that very general magnetic field configurations 
that are inconsequential for the equilibrium state can have profound consequences for local 
WKB perturbations is one of the key conceptual points that emerged in the first accretion 
disk studies of the MRI (Balbus & Hawley 1991, 1992). In this weak field limit, the time 
dependence of the unperturbed magnetic field is a nonissue, which opens a broad range of 
problems to analysis. 

The literature on the hydrodynamical stability of stellar differential rotation is vast. 
The magneto-hydrodynamical stability of differential rotation is a much more specialized 
topic (e.g., Mestel 1999), and has naturally tended to emphasize the direct dynamical forces 
associated with the field itself. The Sun and many other stars are expected to possess 
magnetic fields buried deep inside their radiative zones, as remnants of their complex 
formation history. Their strength is not well known. It is well known, however, that the 
time for buried magnetic fields to diffuse out of the solar interior is very long indeed (see, 
e.g., Parker 1979). Balbus & Hawley (1994) and Balbus (1995; hereafter B95) have studied 
the linear, adiabatic MRI in a stably-stratified stellar system, and have noted that the 
strongly-restoring buoyant forces limit the instability to displacements lying only within 
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spherical shells. The Brunt- Vaisala frequency N in stars is generally orders of magnitude 
larger than the rotation frequency fl (a value of 10'^ is typical of the solar radiative zone), 
whereas is at best comparable to Q in a disk. Even more importantly, the most unstable 
adiabatic MRI modes in a disk are always in the mid-plane (this maximizes the component 
of the displacement along the angular velocity gradient), and these are insensitive to the 
vertical stratification profile. 

Goldreich & Schubert (1967; hereafter GS) showed that the stabilizing effects of entropy 
stratification can be compromised by thermal diffusion. The mechanism is analogous 
to "salt-fingering" in the oceans. In this process, warm salty water overlying cool fresh 
water, naively a stable configuration, is destabilized by heat transfer. Warm fingers of 
salty water, penetrating into the cooler waters below, diffuse heat outward more rapidly 
than they diffuse salt, and thereby loose their buoyancy. In the stellar case, a downwardly 
displaced fluid element is adiabatically heated and is normally warmer than its ambient 
surroundings. This results in a restoring buoyant force. But if thermal diffusion causes 
sufficiently rapid heat leakage, the buoyant force is diminished, and destabilizing angular 
momentum gradients are then able to operate. GS found that not only must the familiar 
Rayleigh stability criterion of increasing angular momentum with increasing axial radius 
be satisfied, the presence of a large thermal conductivity implies that the angular velocity 
must also be constant along cylindrical axes — a far more stringent criterion. 

But great care must therefore be given to apparently small diffusivities when assessing 
the stability of a stellar rotation profile. The GS result holds when the ratio of the thermal 
to viscous diffusivities is sufficiently large. Large compared to what? The answer is not 
unity (Acheson 1978). Rather, it must be large compared to the square of the ratio of 
the Brunt- Vaisala to rotation frequency, a condition not met for the solar radiative zone. 
Acheson (1978) also generalized the study of multi-diffusive modes to non-axisymmetric 
perturbations in a medium with a purely toroidal magnetic field, and found that the 
angular velocity gradient, not the angular momentum gradient, emerged as the rotational 
stability discriminant. This result, it turns out, is very general, extending beyond toroidal 
field geometries (Balbus 1995, 2001). 

In this paper, we solve the problem of multi-diffusive stability to axisymmetric 
perturbations in the presence of rotation, entropy stratification, and a magnetic field of 
arbitrary geometry. We find that the results are sensitive to the triple combination of 
resistivity, viscosity, and thermal conductivity, and that all three must be included in the 
analysis from the very beginning. In what follows, we will refer to such a situation as 
triple-diffusive, whether it is stable or not. 

An outline of the paper is as follows. Section 2 is a derivation of the general dispersion 
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relation described above. Section 3 is a detailed stability analysis. Section 4 applies the 
results specifically to the sun's radiative interior. A more general discussion follows in §5, 
and §6 summarizes our conclusions. 



2. Dispersion Relation 

The MHD equations including the effects of viscosity, resistivity and heat conduction 
take the form 

^^+V.{pv) = 0, (1) 
+ (pv'V)v = - V j^P + ^ j - + (£• B + p (v^v + ^ V(V.t;)) , (2) 



dB 



dt 



V X (v X B) -77VX (V X B), (3) 



1 ^rflnPp- 



■7 



(7 - 1) dt 



^ ^'"1" = XV^T. (4) 



These are, respectively, the continuity equation, the momentum conservation equation, the 
induction equation and the entropy- form of the energy equation (see, e.g., Balbus & Hawley 
1998). Our notation is as follows: v is the flow velocity, p is the mass density, P is the 
pressure, B is the magnetic fleld, $ is the gravitational potential, T is the temperature, 
p is the dynamic viscosity, rj is the resistivity and x is the heat conductivity (which can 
represent thermal or radiative conductivity, depending on the problem at hand). In what 
follows, we write the kinematic viscosity coefficient v — pj p. Bulk viscosity effects are 
neglected. The adiabatic index of the gas, denoted 7, is 5/3 for a monotomic gas with 
neghgible radiation pressure. We have ignored the spatial dependence of the diffusion 
coefficients /i, r/ and X) which is appropriate for a leading order WKB analysis. We have 
also neglected the resistive and viscous dissipation terms in the entropy equation, which are 
also higher order terms. The validity of this approximation is examined in Appendix A. 

We work in cylindrical coordinates {R,(j),Z). We consider axisymmetric Eulerian 
perturbations (denoted by a prefix 6) with WKB space-time dependence exp[i(fe • r — out)], 
where k = {kR,0, kz). The basic state magnetic field, allowed to have any geometry, is 
assumed to be weak enough that it does not affect the basic state configuration, i.e. weak 
compared to both rotation and pressure gradients. The basic state rotation is given by 
f2 = (0, 0, n(P, Z)) along the Z-axis. We neglect in our analysis weak circulations such as 
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those induced by a finite viscosity or any meridional circulation (when considering stellar 
applications) . 

Using the Boussinesq approximation, the leading order WKB terms become after 
linearization 

kRSvR + kzSvz = 0, (5) 



(^—ioi) + uk^^ 8vr 



+ i/k"^^ 5v(j) + 5vr 



H dP - 2ildv^ -— + - — 

X {B^ 5B^ + Bz 5Bz) - ^Bz 5Br = 0, 

47rp 

1 d(R^n) ^ an , ^sb^ 

^ ' +5vzRT:^-ik'B^ = 0, 



R dR 



dZ 



A'Kp 



(6) 
(7) 



(y—iu + vk^^ 5vz 



ikz SP 5p dP ikz 



p p^ dZ ~'~ 47rp 

X {B^ 6B^ + Br 6Br) - 6Bz = 0, 



(^-icu + r]k^^ 5Br - ik • B6vr = 0, 

(-iu + ri^) 6B^ - 6Br - 5Bz R^ - ik - B 5v^ ^ 0, 

(-ia; + Tjk^^ 6Bz - ik • B6vz = 0, 

ia;^^ + (dv . V) InPp-^ = (7 - 1)^^. 
p P p 



(8) 

(9) 
(10) 

(11) 
(12) 



The gravitational potential has dropped from the problem because we are interested in 
wavelengths much shorter than the Jeans wavelength at which self-gravity would become 
important. The relation ST = —T{5p/p) has been used in equation (12) above. 

Solving for the eight 5-unknowns in these eight equations, we obtain the following 
dispersion relation (in compact form) 



— (VP) V In Pp-^ 
IP 



-,v{R^n^] 



An'^{k • VA^i^e = 0,(13) 



where 



VA^B/^finp, k'^^k\ + kl, ujl^^^ujbUJv-{k-VAy, ul ^ ul - {k • va^, 
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ujh = cu + irjk'^ , LUy = LU + iuk'^ , 00^ = 00 + 



kn d d 



kzdZ dR, 



By substituting a = —icu, this 5th order dispersion relation can be written in developed 
form 



where 



Oocr + Oicr + 020" + 0317 + 04(7 + 05 = 0, 



ao — k /k^i 

2uk'^ + 2r]k'^ + ^k'^ 



(14) 

(15) 
(16) 



«2 = T:r 



u^k^ + rj^k^ + 4z/r^r + 2i/^r + 2r]ik'' + 2{k ■ Va) 



IP 



-VP V In Pp- 



(17) 



«3 - T2- 

n7. 



2r]iy'k^ + 2ur]^k^ + u^^k^ + rf^k"" + 4i/?7^r + 2{vk' + r]k' + ^k') {k ■ va) 



-{2r]k'^ + uk'^) 



—VPVlnPp-'^ 

IP 



, (18) 



04 



kz 



2r]iv'e + 2vr]\k'' + r]'y'e + 2(1/77^ + i/^F + r]ik''){k • Va) + {k • v^)' 

1 



-(2i/?7A;^ + 77^A;^ + (fc.t;A)') 



IP 



-VPVlnPp- 



-{2r]^k^ + T]'k^ + {k-VA 



-m\k'VAf, (19) 



«5 = ^ [^^'^'^'° + 2ei^77A;6(fc . vaT + Ck^{k - vaY 



-{urj'e + rjk'{k'VAr) 



—VPVXnPp-'^ 

7P 



-{irfk' + ik\k-VAf) 

-m^{k ' VAf^k"", 



-^V{R'n^) 



(20) 
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where, for conciseness, we have introduced the heat diffusivity 

7-lT 
^ 7 P^' 

The above dispersion relation reduces to several previously established relations in the 
appropriate limits. Taking v — r] = x = ^ (diffusion-free limit), one recovers the result of 
B95 (his Eq. [2.4]), which is a generalization of the work of Balbus & Hawlcy (1991; 1994) 
on the MRI to general rotation laws (i.e. non-constant on cylinders). Taking = f] = Q 
(hydrodynamic limit), one recovers exactly the result of Goldrcich & Schubert (1967; their 
Eq. [32]; see also Fricke 1968). Finally, taking v = 7] = {) but x finite, one recovers the 
result of Balbus (2001; eq. [29]) provided that the isotropic thermal conductivity limit 
is substituted in his dispersion relation. This is effected by setting (fe • 6)^ — > k"^ and 
PlnT — > in the Balbus (2001) dispersion relation. (See also Urpin & Brandenburg 1998). 

3. Stability Analysis 

The complexity of our dispersion relation makes it difficult to derive general necessary 
and sufficient conditions for stabihty. The single-diffusive case (Balbus 2001) is amenable 
to a Routh-Hurwitz (RH) analysis, but we have found this method impractical for the 
multi-diffusive case of interest here. Instead, we derive a series of necessary conditions for 
stability. As we shall see, they are stringent enough to provide very useful limits on the 
maximum allowable level of stable differential rotation. 

To simplify the analysis further, we first consider the inviscid and perfect-conductor 
limits separately. This allows us to reduce the problem to two separate double-diffusive 
situations. As we shall see below, this separation has to be made with some caution. Our 
stability analysis then proceeds as follows. Given that > in Eq. [14], if any of the five 
other tti coefficients is negative, there will be at least one unstable root (i.e. one with a 
strictly positive real part) to the dispersion relation. A necessary condition for stability is 
thus that all the be positive. The requirement ai > is trivially satisfied, so we focus on 
the criteria corresponding to the positivity of the last four coefficients below. Notice that 
all the coefficients 02-05 possess first terms in bracket, oc k^/k^, that are strictly positive. 
These terms represent the systematically stabihzing effects of diffusion processes, {u, rj, 
x), and magnetic tension, {k • Va)^, on very smah scales, where they become dominant. 
Quite generally, we can focus our attention on (at least somewhat) larger scales for 
which the coefficients 02^05 depend on terms combining the effects of differential rotation, 
stratification, diffusion and magnetic tension, and for which positivity is not guaranteed. 
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In addition, we have found that it is possible to carry out a partial RH analysis. 
The first non-trivial determinant required to be positive by the RH criterion is 
det{2) = aia2 — doa^ > 0. This comes as an additional necessary condition for stability. We 
have found that this requirement is important to our conclusions on differential rotation. 



3.1. Perfect-Conductor Limit {rj — > 0) 

This is the case closest to the analysis of Goldreich & Schubert (1967) and Pricke 
(1968), who focused on double-diffusive hydrodynamical stability in the presence of viscosity 
and heat diffusion. 



3.1.1. Requirement a2 > for stability when 7] ^ 

^Prom the structure of the coefficient 02 (involving pure rotational and stratification 
terms), we recognize a stability criterion related to diffusion-free, hydrodynamical modes. 
On large enough scales for the stabilizing effect of the first bracket term in 02 to be 
unimportant, the criterion 02 > becomes 



7P 



-VPVln Pp- 



> 0. 



(21) 



As discussed by B95, this inequality translates into the classical Solberg-H0iland criteria 
(see, e.g., Tassoul 1978) 



dz) 



-- L(VP).VlnPp-- + i^^>0, 
7P R'^ oR 



1 dR^n^ din Pp- 



1 dR'^n^ din Pp- 



R^ OR 



dZ 



i?3 dZ 



dR 



> 0. 



(22) 
(23) 



For a star with spherically-symmetric isocontours of density, p, and pressure, P, these 
criteria take a more familiar form when both cylindrical {R, Z) and spherical (r, 6) 
coordinates are used 



A^^cot^ 



de 



>0, 



(24) 
(25) 
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where the square of the Brunt- Vaisala frequency is defined by 



2 1 dPdinPp-^ 



N' = ^ (26) 

7P or or 

and the square of the epicychc frequency is defined by 

" ~ OR ■ 

Equation (25) shows that to avoid rotational instabihty in a stratified star, the specific 
angular momentum must not decrease away from the rotation axis within a spherical shell. 
While these criteria are necessary and sufficient to guarantee stability in the diffusion-free, 
hydrodynamical limit, they become only necessary conditions for stability in the broader 
context studied here. 



3.1.2. Requirement > for stability when rj ^ 

^From the structure of the coefficient 03, we recognize a stability criterion related 
to hydrodynamical modes influenced on all scales by viscosity and heat diffusion. This 
essentially is the constant term of Goldreich & Schubert's third-order dispersion relation. 
On large enough scales for the stabilizing effect of the first bracket term to be unimportant, 
the criterion 03 > becomes 



—VPVlnPp-^ 
IP 



^,v{R^n^) 



>o. 



(28) 



where the Prandtl number is defined by 



7z/ P 



This condition can be rewritten 



(29) 



where 



fdPdlnPp-^ ^dPdlnPp- 
~p VdZ dR ^ 



R 



1 dP din Pp-^ 
'~p'dR dR ' 



dR dZ 



1 d{R^n^^ 



i?3 dZ 



1 d{R^n^^ 



R^ dR 

1 dPdlnPp-'^ 
~p'dZ dZ ■ 



>0, 



(30) 
(31) 
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As a quadratic in x — kji/kz, this inequality is satisfied if and only if (i) the sum of the 
constant coefficient and the coefficient for x"^ is positive, and (ii) the discriminant of the 
polynomial in x is negative. By using the vorticity relation for the basic state (e.g. Tassoul 
1978) 

dZ p^[dRdZ dZdRj' 



to simplify condition (ii), the two criteria can be written 



+ Nl) + 



'b? or 



>0, 



(33) 



(1 - e.) 



B? dZ ] ^^dZ 



1 d{R^n'^)d In Pp-^ 1 d{R'^n'^) din Pp- 



OR 



dZ 



dZ 



OR 



< 0.(34) 



For a star with spherically-symmetric isocontours of density, p, and pressure, P, these 
criteria take the following form when both cylindrical {R, Z) and spherical (r, 6) coordinates 
are used 

e^N^ + >0, (35) 



(1 - e.) 



dZ 



dd 



< 0. 



(36) 



Note that in the limit — > 0, these two conditions reduce to the result found by Goldreich 
& Schubert (1967): necessary conditions for stability are (i) that the specific angular 
momentum does not decrease with distance from the rotation axis and (ii) that the rotation 
law be constant on cylinders. Achcson (1978) has pointed out, however, that the limit 

^ Q must be carefully taken because in some stars, in particular the Sun, the product 
ei,N'^ is not necessarily small compared to the rotational terms. We discuss this hmitation 
further in §4. Interestingly, for e,^ = 1, the diffusion-free result (§3.1.1) is recovered. 



3.1.3. Requirement 04 > for stability when 77 ^ 

iFiom. the structure of the coefficient 04 (note in particular the — 4r2^(fc • Va)^ term), 
we recognize a stability criterion related to diffusion-free, MHD modes. On large enough 
scales for the stabilizing effect of the first bracket term to be unimportant, the criterion 
04 > becomes 



IP 



-VPVXnPp- 



R^ 



V{R'Q.') + AO.' 



> 0. 



(37) 
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Note the additional assumption made to derive this inequahty: the weak magnetic field 
must be strong enough that magnetic tension forces are important on scales larger than 
those on which dissipation stabilizes all perturbations (e.g. (fc • va)^ ^ rfk^). This still 
leaves a comfortable range of dynamically interesting magnetic field strengths, as is shown 
for the specific Solar case in Appendix A. As described by B95, the above inequality 
translates into 

1 

(VP).VlnPp-^ + — — >0, (38) 

7P a m it 



dn'^dhiPp-^ dQ^ din Pp- 




dZ dR 



> 0, (39) 



i.e. forms similar to the classical Solberg-H0iland criteria but with gradients of angular 
velocity replacing the traditional gradients of specific angular momentum. For a star with 
spherically-symmetric isocontours of density, p, and pressure, P, these criteria take the 
following form when both cylindrical {R, Z) and spherical (r, 6) coordinates are used 

TV^sin^cos^— > 0. (41) 

Note the important new element introduced by the presence of a magnetic field: to 
guarantee axisymmetric stability in a stably-stratified star, the angular velocity must not 
decrease with distance from the rotation axis within a spherical shell. While the two criteria 
above were necessary and sufficient conditions for stability in the diffusion-free study of 
B95, they become only necessary conditions for stability in the broader context studied 
here. 



3.1.4. Requirement 05 > for stability when rj ^ 

/^From the structure of the coefficient 05, we recognize a stability criterion related to 
MHD modes infiuenced by viscosity and heat diffusion, not just on small scales. On large 
enough scales for the stabilizing effect of the first bracket term to be unimportant, the 
criterion > becomes 

> 0. (42) 



R^ 



V{R^n') + 40' 
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Note the remarkable property of this stabihty condition: it is independent of the 
stratification term (which drops out of the dispersion relation in the limit — > 0). (This 
property would not be exactly conserved in a triple-diffusive system.) To be satisfied for 
any combination of ku and kz, it requires the following two conditions 



(43) 



d\nR 



(.§)%o, 

so that marginal stability is possible only for a rotation law that is constant on cylinders 
(dn/dz = 0). 

In a realistic situation, however, even a fully ionized plasma will possess a finite 
resistivity. One should include resistivity in the analysis, no matter how small it is, because 
the stratification term can generally be much larger than the rotational term in a stellar 
context. The situation is intrinsically triple-diffusive because both viscosity and resistivity 
affect the momentum of displaced fluid elements. This makes the analysis of this branch of 
the dispersion relation more complicated, and it appears that a simple stability criterion 
independent of k and va does not exist in general. 



3.1.5. Requirement det{2) > for stability when rj ^ 

Necessary conditions for stability can be made more stringent by the additional 
requirement det{2) = 0102 — ao^s > 0. This is one of the five RH determinants that must 
be positive. 

Keeping all the terms in the coefficients Oq to 03 when calculating det{2), we note 
that several are strictly positive and become negligibly small on large enough scales. The 
condition det{2) > can thus be reduced to 



-(1 + e.) 



—VPVln Pp-^ 
7P 



2e. 



-P(i?^Q^) 



> 0. 



(45) 



Note the different structure of this inequality as compared to Eq. (28). The factor ei, is now 
in front of the rotational term (second bracket). In general, for stars, rotational effects are 
weak compared to the entropy stratification term (first bracket) and <^ 1, so that one is 
tempted to reduce the above expression to the stratification term only. This, however, is 
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incorrect because terms of order N'^/Q'^ emerge, and are important. Using the equivalence 
of Eq. (45) with Eq. (28), if e^, is replaced by (l + e,^)/2e,y in the latter, the analysis proceeds 
as in §3.1.2 and we obtain the following two necessary conditions for stability 



(1 + 6.)(7V| + 7V|) + 



2e^ d{R^n ^^ 
^ dR 



>0, 



(46) 



' 1 djE^n^ 

dZ~ 



+ 



8e,{l + e,) dP 
7P dZ 



1 d{R^n'^)d\nPp- 



1 d{R'^^V)d\YiPp- 



i?3 dR 



dZ 



R3 dZ 



dR 



< 0.(47) 



For a star with spherically-symmetric isocontours of density, p, and pressure, P, these 
criteria take the following form when both cyhndrical (i?, Z) and spherical (r, 9) coordinates 
are used 



(1 - e,y R 



2 I 



dZ 



[1 + €^)N'^ + 2e^K^ > 0, 
- 16e^(l + e^)A^^Qcot^ 



< 0. 



(48) 
(49) 



Note that these conditions differ from those obtained by requiring that 03 > in § 3.1.2. 
The first of the two conditions above requires, in the limit — * 0, that the stratification be 
stable, this time independently of the rotational stability (measured by n^). The second of 
the two conditions above is more stringent than the corresponding one in § 3.1.2 because 
it involves a first term that cannot be negative (or, equivalently, is stabilizing). Only the 
second term can, and its absolute value is a factor of two larger in the above inequalities as 
compared to those we derived in §3.1.2. 

It is worth emphasizing that the combination of requirements 02 > (§ 3.1.1), as > 
(§3.1.2) and det{2) > (above) constitute necessary anrf sufficient stability conditions for 
the double-diffusive hydrodynamical problem considered by Goldreich & Schubert (1967) 
and Pricke (1968). One easily shows that an hydrodynamical system satisfying all these 
conditions also satisfies the RH criterion for the third order dispersion relation of the purely 
hydrodynamical problem. 



3.2. Inviscid Limit (z/ — > 0) 

For the most part, the stability analysis in the limit i/ — > proceeds in a manner 
very similar to the limit 77 — > 0. Unless complications arise, we directly list the necessary 
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conditions for stability in a convenient form. 



3.2.1. Requirement 02 > for stability when u ^ 

Like in the perfect-conductor limit, the necessary conditions for stability are equivalent 
to the classical Solberg-H0iland criteria (see §3.1.1). 



3.2.2. Requirement > for stability when — > 
On large enough scales, the requirement 03 > becomes, in this case. 



—VPVlnPp-'^ 
IP 



- (2e^ + 1) 



>0, 



(50) 



where the "Acheson number^" is 



77/ P 



(51) 



By similarity with the analysis in §3.1.2, we deduce the following necessary conditions for 
stability 



—{Nl + Nl) + — -^r > 0, 

l + 2e„^ ^ ^' R^ dR 



(52) 



(1- 



2e 



V \2 



1 diR'^n') 



l + 2e^' 



7p4 + 2e/9Z 



i?3 dZ J 

" 1 d{R'^Q'^)d In Pp-^ 



1 d{R'^n^)d\nPp- 



R^ dR 



dZ 



i?3 dZ 



OR 



< 0. (53) 



^We were not able to find an existing terminology for this dimensionless number in the literature. Acheson 
(1978) appears to have been the first to recognize the importance of this dimensionless quantity for the 
problem of differential rotation in magnetized and stratified fluids. This number can be expressed, in a 
rather indirect way, as the ratio of the Prandtl number to the magnetic Prandtl number. 
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For a star with spherically-symmetric isocontours of density, p, and pressure, P, these 
criteria take the following form when both cylindrical {R, Z) and spherical (r, 9) coordinates 
are used 



2e. 



2e. 



R 



dQ 
'dZ 



2\ 2 



2e, 



1 + 2e^ 



)N'^ncote 



< 0. 



(54) 
(55) 



In the limit — >^ 0, these two conditions are still consistent with the result of Goldreich 
& Schubert (1967), even though these authors focused on the viscous case, rather than the 
resistive one. The two necessary conditions for stability are (i) that the specific angular 
momentum docs not decrease with distance from the rotation axis and (ii) that the rotation 
law be constant on cylinders. The danger in using the limit — > 0, as noted by Acheson 
(1978) for the viscous case, carries over to the resistive case in the sense that the product 
e^iV^/n^ is not necessarily small compared to unity (see §4). This time, the diffusion- free 
criteria (§3.1.1) are not exactly recovered when = 1. 



3.2.3. Requirement > for stability when — > 

Like in the perfect-conductor limit, the necessary conditions for stability are the 
"modified Solberg-H0iland" criteria derived by B95 (see §3.1.3). 



3.2.4. Requirement 05 > for stability when — > 
On large enough scales, the requirement > becomes 



IP 



-^P(i?"Q2)+4Q2 



> 0. 



(56) 



This condition can be rewritten 




kR 
kz 



'dPdln Pp-^ ^dPdln Pp-^< ' 



dZ dR 



dR dZ 



-R 



'dZ 



dlnR 



> 0. 



(57) 



-16- 



Proceeding as before, the two criteria can be written 



dlnR 



>0, 



(5J 



(1 - e,r R 



dZ 



+ 



dfl"^ dlnPp- 
dlnR dZ~ 



R 



dil^dln Pp- 
'dZ dR~ 



< 0. 



(59) 



For a star with spherically-symmetric isocontours of density, p, and pressure, P, these 
criteria take the following form when both cyhndrical (i?, Z) and spherical (r, 6) coordinates 
are used 



..2 dn^ 



>o. 



R 



'dZ 



on 
'ae 



< 0. 



(60) 
(61) 



In this case, in the hmit — > 0, the first of these conditions differs from the result 
emphasized by Goldreich & Schubert (1967): a necessary condition for stabihty is that 
angular velocity (not specific angular momentum) does not decrease with distance from the 
rotation axis. The second condition, requiring the rotation law to be constant on cylinders 
for marginal stability, remains the same, however. The diffusion-free result (§3.1.3) is 
recovered for e„ = 1. 



3.2.5. Requirement det{2) > for stability when i/ — > 

In this hmit, on large enough scales, the condition det{2) > reduces to the simple 
expression 



7P 



-VPVlii Pp 



-7 



>0, 



(62) 



because the various rotational terms cancel out exactly. 
This condition can be rewritten 
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knll / dP d In Pp-^ dP d In Pp-^ 
[jp [dZ dR ^ dR dZ 



) 



+ Ni>0. 



(63) 



Proceeding as before, the two criteria can be written 



( 



2 



> 0, 



R 



dZ 



<0, 



(64) 



so that marginal stabihty is possible only for a rotation law that is constant on cylinders 
[dQ/dZ = 0). Note that here, the diffusion coefficients 1] and x do not explicitly appear. 
The necessary condition for the rotation law to be constant on cylinders is therefore 
significantly more stringent than the other i/ ^ conditions we have derived thus far. 
Because it relies on the exact cancellation of rotational terms and because it was derived in 
the inviscid limit, however, this result may not strictly hold in a more realistic triple-diffusive 
situation, when the fluid possesses a finite, even if small, viscosity. We revisit this issue 
below when we discuss numerical solutions to our dispersion relation. 



It is significant that in both the inviscid and perfect-conductor double-diffusive limits, 
irrespective of the relative strengths of the various diffusion processes involved (within 
the limitations of our dispersion relation; see Appendix A), rotation must be constant 
on cylinders for stability. In the perfect-conductor limit, this comes from requiring that 
as > 0. In the inviscid limit, this comes from requiring that the second determinant in the 
RH analysis be positive {det{2) > 0). In both cases, we have also recovered as a necessary 
condition for stability the result of Balbus & Hawley (1994) and B95: angular velocity must 
not decrease with distance from the rotation axis within a given spherical shell for stability 
in the stellar context. Consequently, in both double-diffusive limits, a marginally stable 
rotation law must at the same time be constant on cylinders and constant within spherical 
shells. This is achieved only by solid body rotation.^ 

We have cautioned that using one or the other double-diffusive hmit can be somewhat 



^It is important to note here that our focus is on "negative" differential rotation, i.e. differential rotation 
such that dfl/dr < or dQ./d9 < (in terms of the spherical coordinates r and 6). A "positive" differential 
rotation with, for instance, dQ/d9 > and a rotation constant on cylinders is perfectly stable according to 
the criteria we derived. 



3.3. Implications 



-18- 



misleading because the addition of a third, even if small, diffusion process should modify the 
conditions for stability. Because constant rotation on cylinders is related to double-diffusive 
modes of the type described by Goldreich & Schubert (1967) or to corresponding magnetized 
modes in our dispersion relation, we intuitively expect the double-diffusive results to be 
good first approximations as long as focusing on the two largest diffusion processes is 
justified. It should be so when there is a well defined hierarchy of diffusion processes (e.g. 
ly <^ T] ^ i). In general, however, it is possible that some finite amount of stable differential 
rotation persists in a fully triple-diffusive situation. We have not been able to establish 
necessary conditions for stability in the general, triple-diffusive case and we have addressed 
this issue by numerically solving the full dispersion relation. 



4. Numerical Solutions for the Sun's Radiative Zone 



The search for numerical solutions to our dispersion relation described in this section 
has three goals: (1) to confirm the results of our double-diffusive analysis, (2) to explore 
the potentially stabilizing role of a third, weaker diffusion process on an otherwise unstable 
double-diffusive situation, and (3) to apply our results to the Sun's radiative zone. 

We adopt a standard model for the current Sun (e.g. Bahcall, Pinsonneault & Basu 
2001; Demarque & Guenther 1991). We focus on the part of the radiative zone, from 
r ~ 0.3 — O.7i?0, in which composition gradients are small. We estimate the values of 
various microscopic parameters relevant to the stability problem. Following Spitzer (1962), 
the (ion-dominated) dynamic viscosity for a hydrogen-dominated plasma is 



rp5/2 

= pi/ ~ 2.2 X 10"^^^ g cm-^ s-\ (65) 

where In A ~ 4 is an appropriate value of the Coulomb logarithm for the Solar interior. The 
resistivity for a hydrogen-dominated plasma is 



77 ~ 5.2 X 10^^^ cm^ s^K (66) 

Radiative heat diffusion dominates over thermal heat diffusion in the solar interior, 
with a radiative conductivity given by (e.g., Schwarzschild 1958) 



_ 16TV 
3kp 



(67) 
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where a is the Stefan-Boltzmann constant and k is the radiative opacity. The corresponding 
radiative diffTisivity, which can be directly compared to the kinematic viscosity and 
resistivity, is given by 



irad— — ^ — -pXrad- (68) 

We hst in Table 1 the values of the density, p, temperature, T, Rosseland-mean opacity, 
K (obtained from a standard opacity table) and all three diffusivities over the region of 
interest in the solar radiative zone. In addition, the values of the Prandtl and Acheson 
numbers, and e^, which are directly relevant to the stability analysis, are listed. Clearly, 
resistive diffusion dominates over viscous diffusion in the bulk of the solar radiative zone, 
by a factor ~ 20 — 30. Note that the radiative kinematic viscosity 



^ 15^^' ^''^ 

where c is the speed of light, makes only a small contribution to the total kinematic viscosity 
in that region (e.g. Goldreich & Schubert 1967). In our specific numerical applications, 
we focus on the region below the convection zone at r ^ 0.7, for which helioseismological 
measurements of a near solid-body rotation are most reliable, with an angular velocity 
Q — 2.7 X 10~^ rad (see, e.g., Charbonneau, Dikpati & Oilman 1999). We adopt the 
value — 1.3 X 10~^ Hz^ for the square of the Brunt- Vaisala frequency (Demarque & 
Guenther 1991; their Fig. 23), so that the ratio Q'^/N^ = 5.6 x 10"^ in our models. 

It is important to note at this point that, in the Sun's radiative zone, the Prandtl 
number, e^, is comparable to D,^/N^ while the Acheson number, e^, is systematically 
» Q?' /N"^. Indeed, Acheson (1978) has pointed out that, since the double-diffusive stabihty 
criteria of Goldreich & Schubert (1967) are typically of the form e^N'^ + > (see §3.1.2), 
substantial differential rotation, for instance a Keplerian-like profile with = — can 
remain stable depending on the exact values of e,^, Vt^ and N"^. This objection could be a 
practical limitation of the results of Goldreich & Schubert (1967; see also Fricke 1968), who 
focused on the idealized limit e^, — >■ 0. Acheson"s point is made even stronger by observing 
that the relevant double-diffusive limit for the Sun is the inviscid one (since X ^ ^ > ^) 
and that » ft^/N^ (see Table 1), which appears to be strongly stabihzing by Acheson's 
arguments. Our extensive stability analysis in the previous section shows, however, that 
in both the inviscid and perfect-conductor double-diffusive limits, any level of negative 
differential rotation is destabilized by a combination of diffusion-free and double-diffusive 
modes, irrespective of the relative strength of viscous, resistive and heat diffusion. This 
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result effectively invalidates Acheson's objection but it remains to be seen what the effects 
on stability of adding a third, weaker diffusion process may be. 

We have numerically solved the complete dispersion relation (Eq. [14]) with parameters 
appropriate for the Sun's upper radiative zone, using the Laguerre algorithm described by 
Press et al. (1992). For definiteness, we have considered only conditions appropriate to the 
specific radius r ~ 0.7 Rq: z/© = 23.6 cm^ s~^, rjQ = 596 cm^ s~^ and Crade = 1-2 x 10'^ 
cm^ s^^. In general, we considered a range of values for the polar angle, 6, in the interval 
[0,7r/2] (pole to equator). We have found it useful to rewrite both the rotational and 
stratification terms appearing in the coefficients 02-05 of the dispersion relation in terms 
of spherical coordinates {r,6,4>). For a spherical star, the stratification term becomes only 
function of N^, 9 and the radial and angular wavevectors, kr and ke. The rotational term, on 
the other hand, explicitly depends on the amount of differential rotation within and between 
spherical shells, that we express as d In Q/ 89 and dlnfl/dlnr, respectively. As we have 
mentioned earlier, we are mostly interested in negative differential rotation, i.e. cases where 
d In fl/ 89 < and/or dlnQ/d\nr < 0. We will consider these two cases separately because 
we expect differential rotation within spherical shells to be destabilized by diffusion-free 
modes and differential rotation between shells to be destabilized by multi-diffusive modes. 

When searching for unstable modes, we vary the wavevectors kr and ke independently 
in the range 



27r ^ 27r 

^ lO-i^i?©' 

The first (large-scale) limit guarantees that we are looking at scales significantly smaller 
than the pressure scale height {Hq ~ O.IRq), while the second (small scale) limit guarantees 
that we are looking at scales significantly in excess of the mean free path for the conditions 
of interest. We have also performed focused searches on small regions of the wavevector 
space corresponding to nearly cylindrical-radial {\kji/kz\ — >■ 0) displacements or nearly 
vertical {\kji/kz\ — ^ 00) displacements. Independently of the values of the wavenumbers, we 
have varied the Alfven speeds vat and va9 in the range 10~^ - 10~'^^RqQ and we have also 
explored cases with vat — and vao — 0. This large range of values for the Alfven speeds, 
independently of the values for kr and ke, effectively provides a search both in magnetic 
field strength and geometry. With this extensive search, we have explored regimes in which 
|fe • t;^! is successively <C and comparable to each of the three dissipation terms (z^/c^, 
rjk^ and Cradk"^)- 



In the inviscid limit {u — 0), we were able to find unstable modes down to values of 
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differential rotation between shells as low as dhiQ/dlnr ^ —0.1 to —0.01, depending on 
the value of the polar angle, 9. In the perfect-conductor limit (77 = 0), we were able to 
find unstable modes down to values as low as d\nQ/d\nr < —0.01 for essentially all polar 
angles. All the unstable modes were of the direct type {Im{a) = 0), as opposed to the 
overstable type {\I'm{a)\ > Re{a) > 0), as expected in the situation of interest, with a 
strongly stabilizing thermal stratification and a fast rate of heat diffusion. That we found 
unstable modes down to such very low values of negative differential rotation between shells 
in both double-diffusive limits is consistent with the conclusions of our stability analysis in 
§3. 

We then turned to the fully triple-diffusive situation and first investigated differential 
rotation within spherical shells (i.e. dlnQ/d9 < only). We expect this type of differential 
rotation to be destabilized by diffusion-free modes, which should be easily identified even in 
the triple-diffusive case. Indeed, we have found that unstable modes were easily identified 
down to levels dlnQ/d9 ^ —0.1 to —0.01, depending on the value of the polar angle, 9. 

We then searched for unstable modes in the presence of differential rotation between 
shells (i.e. dlnQ/dlnr < only). We have been able to identify unstable modes down 
to levels of differential rotation corresponding to dlnQ/dlnr ~ —1.2 to —1.5, depending 
on the value of the polar angle, 9. It became clear during the extensive search required 
to identify these triple-diffusive modes that they are located in a much smaller region of 
the parameter space than the unstable modes we previously identified in the corresponding 
double-diffusive limit. This is as expected if the addition of a third, weaker diffusion 
process stabilizes an otherwise unstable double- diffusive situation. We have confirmed this 
stabilizing effect explicitly by observing that an initially weakly unstable triple-diffusive 
situation slowly makes a transition to stability as the value of the third, weakest diffusion 
coefficient is increased. 

In that respect, it is worth noting that the triple-diffusive situation in the Sun is such 
that the weakest diffusion process, viscosity, is "only" a factor 20-30 times smaller than 
the second weakest one, resistivity, and this may be a source of stabilization with respect 
to small levels of differential rotation between shells. During our numerical exploration 
of the triple-diffusive situation in the Sun's upper radiative zone, we have also noticed a 
rather strong sensitivity of the stability results to small variations (say, x 2) in the values 
of parameters such as the viscosity, i^q, or the Brunt- Vaisala frequency, A^. A complete 
stability analysis for the current Sun is beyond the scope of the present study. It is 
encouraging, nonetheless, that we have been able to identify unstable modes in the presence 
of moderate levels of differential rotation for the current upper radiative zone conditions, as 
it suggests that these modes may have played an important role in establishing the current 
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rotation profile in the Sun. 

5. The Emergence of Solid Body Rotation 

It lias long been realized that the timescale for microscopic viscosity to reduce significant 
levels of differential rotation in the Sun is excessively long {Rq/uq ~ 0.5-1 x 10^^ years; 
see, e.g., Goldreich & Schubert 1967). On the other hand, differential rotation (with 
dlnQ/dlnr < 0) is expected to have been present in the early Sun, because of the likely 
fast initial rotation and of the magnetic spin-down torque externally exerted via the solar 
wind (see, e.g., Sofia et al. 1991 for a review). Consequently, the solid body rotation 
inferred for the Sun's upper radiative zone from seismology (e.g., Kosovichev et al. 1997; 
Schou et al. 1998; Charbonneau et al. 1999) requires a mechanism capable of reducing 
differential rotation much more efficiently. 

Several mechanisms have been proposed and discussed at length in the literature 
(see, e.g., Schatzman 1991 for a review). They include the Goldreich-Schubert-Fricke 
(GSF) instability, the "secular hydrodynamical shear instability" (Zahn 1974) and angular 
momentum transport by internal gravity waves (Press 1981; Kumar & Quataert 1997; 
Talon Sz Zahn 1998; Talon, Kumar & Zahn 2002). None of these mechanisms may be able 
to provide a satisfactory solution, however. The GSF instability alone is not expected to 
bring a system to a state of solid body rotation (but only of rotation constant on cylinders; 
e.g., Goldreich & Schubert 1967). The existence of a secular shear instability has not been 
rigorously proven, but only inferred from heuristic arguments (Zahn 1974; Schatzman 1991), 
and it would operate in the current Sun's radiative zone only in the presence of strong 
shear (Zahn 1993). Finally, while internal gravity waves may transport angular momentum 
efficiently, they are not generally expected to be efficient at mixing elements. The large 
amount of Li depletion at the surface of the Sun and other stars (see, e.g., Ghaboyer, 
Demarque & Pinsonneault 1995a,b) is best interpreted as resulting from turbulent mixing 
in stellar radiative zones, thus favoring the action of instabilities rather than waves to 
explain both mixing and rotational evolution. 

According to our analysis, it is the combination of weak magnetic fields and multi- 
diffusive modes that may render even a small level of (negative) differential rotation 
unstable. Independently of this issue of stability, there are two important aspects of the 
problem that a linear analysis cannot address. First, it is a priori unclear whether the 
turbulence resulting from the non-linear development of the various modes described by 
our dispersion relation will drive the system towards a state of marginal stability (i.e. 
near solid body rotation). Intuitively, because the main force balance in a star does not 
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involve rotation (contrary to an accretion disk), we would expect the turbulence to be 
relatively "free" to bring a star close to a state of marginal stability, but this remains to 
be proven. Second, it is also unclear with what efficiency (i.e. on what timescale) that 
turbulence would be able to transport angular momentum and affect an unstable rotation 
profile. Addressing these two non-linear aspects of the problem reliably would require fully 
turbulent numerical simulations. In that respect, it is interesting to note that a preliminary 
two-dimensional investigation by Korycansky (1991) of hydro dynamical double-diffusive 
modes in a specific "equatorial" geometry does indicate that angular momentum transport 
drives the system towards marginal stability. An additional motivation for carrying out 
such detailed numerical simulations of turbulence driven by multi-diffusive, magnetized and 
unmagnetized modes would be to estimate the efficiency of turbulent transport of elements. 

One important physical element that has been neglected in our analysis is the 
stabihzing effect of composition gradients. Goldreich & Schubert (1967) have described, 
in the single-diffusive hydrodynamical limit, how even moderate gradients of chemical 
composition can stabilize significant levels of radial differential rotation. As these authors 
have noted, this effect would be important at radii r ^ O.Si?© in the Sun's core, a region 
that we have excluded from our analysis. While a complete derivation of our dispersion 
relation including the effects of composition gradients is beyond the scope of the present 
study, by analogy with the results of Goldreich & Schubert (1967), one may expect 
significantly stronger levels of radial differential rotation to be maintained in the solar core, 
relative to the rest of the radiative zone. On the other hand, it is possible that the early 
differential rotation in the region currently encompassing the solar core was reduced before 
significant hydrogen burning took place, if multi-diffusive modes were rather efficient early 
on at redistributing angular momentum. This illustrates how the problem of differential 
rotation is intimately linked to that of stellar evolution. 



6. Conclusion 

We have studied the local axisymmetric triple-diffusive stabihty of stratified, weakly- 
magnetized, differentially-rotating fluids. We have established that, in an inviscid or 
a perfectly-conducting fluid, differential rotation is destabilized by a combination of 
diffusion-free and double-diffusive modes, unless rotation is constant on cylinders and 
angular velocity does not decrease away from the rotation axis. We have stressed the 
important role of weak magnetic flelds in establishing these results. We have found that, in 
a more realistic triple-diffusive situation, the weakest diffusion process can sometimes play a 
stabilizing role. While our analysis is rather general, we have discussed a specific numerical 
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application to the Sun's upper radiative zone, which is seismologically known to be rotating 
near solid body rotation. We have found that moderate to strong levels of differential 
rotation, if present, would indeed be destabilized, thus suggesting that magnetized and 
multi-diffusive modes may have played an important role in establishing the current solar 
internal rotation. 
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Appendix A: Limitations 

Al: Validity of the dispersion relation 

We neglected the resistive and viscous dissipation terms in the entropy equation 
(Eq. [4]; see, e.g., Tassoul 1978 or Balbus & Hawley 1998 for complete formulations) and its 
subsequent linearized form (Eq. [12]). We explore the range of validity of this approximation 
here. 

Let us first compare the magnitudes of the three terms we kept in equation (12) 
for entropy perturbations. The ratio of the first term on the LHS to the perturbed heat 
diffusion term is of order 



R,^ "^^/^ ^ I (70) 

' xk''5T/P xk^T/P {i/v)\csk^ {^/u){kH){kXy ^ ' 

where we have used the definition of the heat diffusivity, ^, and we have equated the 
kinematic viscosity coefficient, z/, to the product of the mean free path. A, and sound 
speed, Cg. Note that the perturbation frequency, u, is typically (though not exclusively) 
~ N (the Brunt- Vaisala frequency) in a pressure-supported system and ~ Q (the rotation 
frequency) in a rotation-supported system. In both cases, uj ~ Cg/H, where H is the 
pressure scale-height of the system. Irrespective of the ratio of diffusivities {^/v), our local 
analysis making use of the MHD equations is vahd only for scales much smaller than the 
system's scale-height {kH ^ 1) and much larger than the mean free path {k\ <^ 1). In 
general, the ratio Ri can therefore be ^ 1 or ^ 1 and one must keep the heat diffusion 
term in the perturbed entropy equation to leading order. 
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The ratio of the same first term on the LHS of equation (12) to the (neglected) 
perturbed viscous heating term is of order 

R, ^ ^^^/^ , ^ (71) 

where we have used the relations Sv/cg ~ Sp/p ~ ST/T and \dQ/dlnr\ ~ Q. Since co/Q is 
typically ~ 1 (rotation-supported system) or ^ 1 (pressure-supported system), and kX <^ 1, 
R2 is ^ 1 and neglecting the perturbed viscous heating term is justified to leading order. 

Similarly, the ratio of the first term on the LHS of equation (12) to the (neglected) 
perturbed resistive heating term is of order 



^ u;6T/T u;6T/T 1 

' 6[{7]/47r)\VxBf]/P^ {v/^7r){vA/cs)^k6cs/{csH) {r]/Mv){vAlcsf{k\y ^ ' 

where we have made the additional assumption that the ratio va/cs is locally constant 
in the basic state configuration. This amounts to requiring that the basic state magnetic 
field does not possess strong gradients on scales smaller than the scale height, if, and is a 
reasonable approximation unless one is interested in rather singular basic state magnetic 
field configurations. Irrespective of the ratio of diffusivities, rj/v, our weak field assumption 
implies that Va/cs <^ 1, so that R-^ is also ^ 1 in general and neglecting the perturbed 
resistive heating term is justified to leading order. 

Physically, this hierarchy of heating terms in the perturbed entropy equation can 
be understood by noting that heat diffusion is naturally the most efficient form of heat 
transport, unless the viscosity and/or resistivity coefficients are very much larger than the 
heat diffusivity coefficient (i.e. v ox r] ^ ^). That our assumptions for the leading-order 
perturbed equations remain valid even if i/ or 77 is larger than ^, but not so much as to 
become more efficient at transporting heat than heat diffusion itself, indicates that our 
dispersion relation is able to describe overstable modes in a strongly stratified medium 
where such a hierarchy of diffusivities occurs. 



A2: How weak can weak magnetic fields be? 



In the coefficients 04 and 05 of our dispersion relation, pre-factors multiplying the 
rotational and stratification terms involve sums of dissipation and magnetic tension terms. 
Our double-diffusive stability analysis in §3 assumes that the magnetic field strength, while 
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weak {va <S Cg and va ^ (HRY^'^Q), is strong enough for the magnetic tension term to 
dominate over the dissipation terms in all these pre-factors. We determine the range of 
validity of this assumption here, using numerical values appropriate for the Sun. 

Let ^ denote the largest diffusivity coefficient (in units of cm^ s~^) in the medium of 
interest. For the assumption made in our double-diffusive stability analysis to be incorrect, 
we need k'^v\ <^ $,^k^ for all possible values of k relevant to the local analysis. Since the 
minimum acceptable value of A; is ~ 2tt/Hq (where Hq is the scale height), our assumption 
breaks down for all relevant scales if v\ ^ An'^^'^ / Hq. Using typical values p ~ 1 g cm~^, 
^ ~ 10^ cm^ s~^ and Hq/Rq ~ 0.1 for the Sun's upper radiative zone, this translates 
into a limit on the field strength i? -C 8 x 10^^ G, a very small value indeed. Since the 
above-mentioned pre-factors in 04 and 05 actually involve products of several diffusivity 
coefficients (not just the largest one), the hmit on the field strength will be smaller 
than we estimated by several extra orders of magnitude for the conditions in the Sun's 
upper radiative zone. On the other hand, our assumption of weak magnetic field leads to 
v\ = B^/Anp < {HqRqY/'^Q. < or S < 10^ G. This still leaves a comfortable range of 
field strengths for which the assumptions made in our double-diffusive stability analysis are 
valid for the solar interior. 



REFERENCES 

Acheson, D. J. 1978, Phil. Trans. Roy. Soc. London, 289A, 459 
Armitage, P. J. 1998, ApJ, 501, L189 

Bahcall, J. N., Pinsonneault, M. H. & Basu, S. 2001, ApJ, 555, 990 

Balbus, S. A. & Hawley, J. F. 1991, ApJ, 376, 214 

Balbus, S. A. & Hawley, J. F. 1994, MNRAS, 266, 769 

Balbus, S. A., & Hawley, J. F. 1998, Rev. Mod. Phys., 70, 1 

Balbus, S. A. 1995, ApJ, 453, 380 (B95) 

Balbus, S. A. 2001, ApJ, 562, 909 

Balbus, S. A. 2003, ARA&A, in press (preprint) 

Blaes, O. 2003, in "Accretion Disks, Jets, and High Energy Phenomena in Astrophysics", 
Proceedings of Session LXXVIII of Les Houches Summer School, August 2002, eds. 
F. Menard, G. Pelletier, G. Henri, V. Beskin, and J. Dalibard (EDP Science: Paris 
and Springer: Berhn), astro-ph/0211368 

Chaboyer, B., Demarque, P. & Pinsonneault, M. H. 1995a, ApJ, 441, 865 



-27- 



Chaboyer, B., Demarque, P. & Pinsonneault, M. H. 1995b, ApJ, 441, 876 
Charbonneau, P., Dikpati, M. & Gilman, P. A. 1999, ApJ, 526, 523 

Demarque, P. & Guenther, D. B. 1991, in Solar Interior and Atmosphere, eds. A.N. Cox, 
W.C. Livingston and M.S. Matthews (Tucson: University of Arizona Press), p. 1186 

Fricke, K. 1968, Z. Ap., 68, 137 

Goldreich, P. & Schubert, G. 1967, ApJ, 150, 571 

Hawley, J. F. 2001, ApJ, 554, 534 

Hawley, J. F., Gammie, C. F. & Balbus, S. A. 1995, ApJ, 440, 742 
Korycansky, D. G. 1991, ApJ, 381, 515 
Kosovichev, A. G. et aL 1997, Sol Phys., 170, 43 
Kumar, P. & Quataert, E. J. 1997, ApJ, 479, L51 

Parker, E. N. 1979, Cosmical Magnetic Fields (New York: Oxford University Press) 

Press, W. H., Teukolsky, S. A., Vetterhng, W.T. & Flannery, B. P. 1992, Numerical Recipes 
in Fortran. The art of scientific computing (Cambridge: Cambridge University 

Press) 

Schou, J. et al. 1998, ApJ, 505, 390 

Schwarzschild, M. 1958, Structure and Evolution of the Stars (New York: Dover) 

Sofia, S., Kawaler, S., Larson, R. & Pinsonneault, M. 1991, in Solar Interior and Atmosphere, 
eds. A.N. Cox, W.C. Livingston and M.S. Matthews (Tucson: University of Arizona 
Press), p. 140 

Schatzman, E. 1991, in Solar Interior and Atmosphere, eds. A.N. Cox, W.C. Livingston 
and M.S. Matthews (Tucson: University of Arizona Press), p. 192 

Spiegel, E. A. & Zahn, J. P. 1992, A&A, 265, 115 

Spitzer, L. 1962, Physics of Fully Ionized Gases (New York: Wiley Interscience) 

Talon, S. & Zahn, J.-P. 1998, A&A, 329, 315 

Talon, S., Kumar, P. & Zahn, J.-P. 2002, ApJ, 574, L175 

Tassoul, J.-L. 1978, Theory of Rotating Stars (Princeton: Princeton University Press) 
Urpin, V. & Brandenburg, A. 1998, MNRAS, 294, 399 

Zahn, J.-P. 1974, in Stellar Instability and Evolution, Proc. lAU Symp. 129, eds. P. 
Ledoux, A. Noels and A.W. Rodgers (Dordrecht: D. Reidel), p. 185 



-28- 

Zahn, J.-P. 1993, in Les Houches School (Session XLVII): Astrophysical Fluid Dynamics, 
eds. J.-P. Zahn & J. Zinn- Justin (Elsevier), p. 561 



This preprint was prepared with the AAS IM^jX macros v4.0. 



-29- 



Table 1: Diffusive Conditions in the Solar Radiative Zone 



Radius 


P 


T 


K 


V 


Vj. 






El/ 




{Re) 


(g cm-3) 


(106 K) 


(cm^ g-i) 


(cm^ s~^) 


(cm^ s~^) 


(cm^ s~-^) 


(cm^ S~^) 






r ~ 0.7 


0.2 


2.3 


18 


21 


2.6 


596 


1.2 X 10^ 


2 X 10-6 


5 X 10-5 


r ~ 0.3 


8.5 


6.2 


3 


6 


0.5 


135 


8.1 X 10^ 


8 X 10-6 


1.7 X 10-^ 



